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I. INTRODUCTION 

The goal of this project is to design single-chip low-pass filters with constant group delay in the 
pass band and 60 dB minimum attenuation in the stop band. The desired 3dB frequencies are (in 
Hz): 

2.5, 5, 10, 20, 40, 80, 160, 320 

A filter class that satisfies the constant delay (linear phase) requirement while providing quite a 
narrow transition band is the Bessel-Chebyshev filters [1]. It was found that the 7th order 
Bessel-Chebyshev response satisfied the requirements of the filtCT. The filter transfer function is 


Ha{s) = K 


(52+8.329^)G2+io.572)(^2+i9.3i2) 

s7+2856+37855+3150s4+17325s3+6237052+135135s+135135 


( 1 ) 


where /T is a normalization constant so that //(O) = 1. The 3dB frequency (normalized) is at 12 = 
2.286. The denominator is the 7th order Bessel polynomial. The magnitude and phase plots of 

HJJCi) are shown in figure 1. 



Figure 1. Loss and group delay of the 7th order Bessel-Chebyshev filter 
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The implementation was as a switched-capacitor (SC) network. These are sampled data systems, 
but the signal amplitude is not quantized [2]. The SC technique allows high order filters to be 
integrated on silicon with very good accuracy. The shape of the transfer function depends on the 
capacitor ratios which are quite precisely controlled. The frequency axis can be scaled by changing 
the system clock frequency. 

II. CASCADE SYNTHESIS 

The transfer function in equation (1) can be expressed as the product of lower order functions of s, 
namely bilinear and biquadratic functions. The desired transfer function is obtained when circuits 
realizing these simpler functions arc cascaded. The rest of this section describes the design steps. 

A. Factorization of the transfer function: 

The numerator of (1) is already factorized, the denominator can be factorized using the 
Newton-Raphson root- finding algorithm. The poles and zeros of (1) are given below; 


Po= -4.9717869 
Pi = -2.6856769-1- j5.4206941 
-4.0701392-h j3.5171760 
P 3 = -4.7582905-h j 1.7392861 

Zq= oo 

zi= J8.329 
22 = j 10.57 
23 = jl9.31 


P 4 = -2.6856769- j5.4206941 

P 5 = -4.0701392- j3.5171760 

P 6 = -4.7582905- j 1.7392861 

24 = -j8.329 
25 = -j 10.57 
26 = -jT9.31 
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Therefore, the transfer function can be decomposed into one bilinear and three biquadratic 
functions. 


B. Pairine of poles and zeros: 

After we pair z, and p, to form the bilinear factor, there remains the issue of pairing the remaining 
three complex conjugate pairs of poles with the three complex conjugate pairs of zeros. 

Out of the many possible pairings of poles and zeros some result in transfer functions that are less 
sensitive to element nonidealities of the actual circuit implementation. It is known that the 
magnitude of a bounded transfer function is not sensitive to element variations at the frequencies 
where it achieves its maximum [3]. Motivated by this property, we can define our objective in the 
pairing of poles and zeros as that the magnitude of the transfer function of each of the building 
blocks must be as flat as possible in the pass band of the filter [4]. 

The flattest blocks occur when the pole-zero pairs are (pj, z{), (P2, Z2), (P3, Z3). The Q factors of 
each block (pole) are as follows: 

01 = 1.126, (pi) 

02 = 0.661, (P 2 ) 

03 = 0.532, (P 3 ) 

C. Ordering of the blocks: 

As it is with the pairing of poles and zeros, different orderings of the blocks result in different 
sensitivities to the nonidealities in the implementation. The objective is again to have the flattest 
possible response in the pass band at the output of the each block . It turns out that the best 

ordering is (po> Zq)— >(Pi, zi)-»(p 2 , Z3). Let us define the transfer functions for each 

block: 
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( 2 ) 


Hoa(s) = ko 
Hia(.s) = ki 

H2a(s) = k2 

H3a(.s) = ky 


1 

^+4.9717869 

(s2+8.3292) 


(52+5.3713538^+36.5967852) 

(52+10.57^) 

( 52 + 8 . 14027845+28.9365601) 

(52+19.312) 

( 52 + 9 . 5 1 659065+25.6664902) 


(3) 

(4) 

(5) 


D. Mapping to the discrete time domain: 

We have dealt with the continuous-time frequency transfer function until now. We use the bilinear 
mapping to the discrete-time complex frequency, z, domain: 


5 = ljtl 
^ 7z+l 


( 6 ) 


where T is the sampling period and the continuous- and discrete-time complex frequencies are 
respectively, 


From (6) we obtain 


s = 2: + j £2 

z = reJ^ 


coT/2 = arctan(f2T / 2) 


(7) 


If or is very small then the relationship between co and £2 is almost linear. 


Applying (6) to (2) through (5) gives the transfer functions to be realized. The numerical values 
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depend on the particular value of T. 


E. Filter sections: 

We realized the filtCT as a cascade of one bilinear and three biquadratic sections. 


The circuit diagram of the bilinear filter section is shown in figure 2. 



Figure 2. The bilinear SC section. 

The transfer function in terms of the capacitor values is given below (C^ is set to 1): 


(1+C2)z-1 


( 8 ) 


The circuit diagram of the high-Q biquadratric filter section is shown in figure 3. 
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The circuit diagram of the low-Q biquadratric filter sections is shown in figure 4. 



Figure 4. The low-0 biquadratic SC section. 
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The transfer function in terms of the capacitor values is given below (Q and Cg are set to 1): 


( 2 ) = . C5Z^+(CiC3+-2Cs)z+C5 

’ (1+C4)z2+(C2C3-C4-2)z+1 


( 10 ) 


The calculated capacitor values are such that when the clock frequency (1/T ) is 10 kHz, the 3 dB 
frquency is 40 Hz. The capacitors are also scaled for maximum dynamic range and minimum size. 
The resulting values are shown in Table 1. 


TABLE 1 

Capacitor values for 3 dB frequency at 40 Hz with 10 kHz clock 



Capacitor 

Value (in unit capacitors) 

Bilinear Section 

Q 

35.58923 


Cl 

1.0 


C2 

2.0 


^3 

2.0 

High-0 Biquadratic Section 

Ca 

7.28401 



31.98662 


Cl 

1.0 


^2 

1.0 


C3 

1.0 


Q 

13.34987 



16.40688 

1st Low-Q Biquadratic Section 

Ca 

25.58748 


Cb 

10.68347 


Cl 

1.0 


C2 

1.0 


C3 

1.0 


C4 

1.0 


C5 

2.90374 

2nd Low-Q Biquadratic Section 

Q 

33.72466 



13.62551 


Cl 

1.0 


Cl 

1.0 


C3 

1.50305 


C4 

1.50305 


C5 

1.0 
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F. Filter sections with reduced capacitor spread: 

A second filter was designed to realize much lower passband edge fi^uencies with moderate clock 
speeds. We used capacitive voltage dividers to reduce the capacitor spread to realizable levels. The 
filter has 3 dB loss at 2.5 Hz when operated with 10 kHz clock. 

The circuit diagrams for the filter sections with capacitive voltage dividers (T-sections) are shown 
in Figures 5 through 7. 
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Figure 5. Bilinear section with voltage dividers. 
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Figure 7. Low-Q biquadratic section with voltage dividers. 


The capacitor values are given in Table 2. 
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TABLE 2 

Capacitor values for 3 dB frequency at 2.5 Hz with 10 kHz clock 


Value 

Capacitor with T 


Bilinear Section Q 23.19570 

Cl 1.0 

Ci5 1.0 

Cjr 23.19570 

C 2 1.43863 

C 25 1.43863 

C 2 T 23.19570 

C 3 1.43863 

C 35 1.43863 

C-ij 23.19570 

High-j2 Biquadratic Section C^ 4.87619 

Cb 21.51080 

Cl 1.0 

C\s 1.0 

CiT 21.51080 

C 2 1.0 

C 25 1.0 

C 2 T 21.51080 

C3 1.0 

C35 1 .0 

C 2 T 21.51080 

C 4 9.08509 

C5 11.32699 

1st Low-g Biquadratic Section C^ 19.25835 

Cb 8.38634 

Cl 1.0 

Cis 1.0 

Cir 19.25835 

C 2 1.02497 

C 25 1.02497 

C 2 T 19.25835 

C 3 1.0 

C 35 12.0 

C 37 19.25835 

C 4 1.02497 

C 5 2.17828 


Value 
without T 


584.42578 

1.0 


2.0 


2.0 


114.64308 

505.73633 

1.0 


1.0 


1.0 


213.59784 

266.30664 

409.40088 

178.27975 

1.0 


1.0 


1.0 


1.0 

46.30406 
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2nd Low>^ Biquadratic Section 

Ca 

26.28833 

539.59912 



14.85170 

152.42438 


Cl 

1.18146 

1.0 


Ci 5 

1.18146 



Cyr 

26.28833 



Cl 

1.0 

1.0 


C2S 

1.5 



Cir 

26.28833 



C3 

1.70085 

1.0 


Css 

1.70085 



Cyr 

26.28833 



C4 

2.0 

1.0 


C5 

1.02569 

10.52675 


As can be seen from Table 2, the overall capacitor area is reduced by a power of 2. 


III. FABRICATION AND TEST RESULTS 

Both of the designs were fabricated as a single chip using 5-|xm double-poly CMOS process 
provided by MOSIS, a federal government supproted agency. We obtained the samples in March 
1988. Table 3 gives some of the physical properties of the fabricated chips: 


TABLE 3 

Some Physical Properties of the Fabricated Filter 


Technology: 

p-well CMOS, double poly 

Minimum feature size: 

5 |J.m 

Chip size: 

6800 |im X 4595 pm 


(265 mil X 179 mil) 

Unit capacitor: 

0.978 pF 

Power supplies: 

±5V 

Power consumption: 

280 mW 

Clock input: 

±5 V square-wave with 50% duty cycle 

Package: 

40 pin ceramic DIP 


The chip phograph is shown in Figure 8. 


The yield was somewhat low, around 50%. This could have been due to the large size of the chip 
by typical MOSIS standards. The frequency domain measurments were done using a network 
analyzer. The magnitude and phase response plots are provided at the end of the section. The 
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results for both sets of filters are in close agreement with the design values. Due to the limitations 
of the network analyzer the plots do not show the response in the 0 to 5 Hz range. Because of this 
limitation, the plots for the filter with the voltage dividers were obtained with 160 kHz clock 
frequency. 

The manual testing at DC showed that the filter with voltage dividers had quite high output offset 
voltage, around 0.8 V. But the other filter had much lower output offset, around 20 mV. We have 
found that the offset voltage can be reduced by adjusting the rising and falling times of the input 
clock signal. 

The filter that was designed to give 40 Hz passband edge at 10 kHz clock operated without 
problems for clock frequencies in the range 625 Hz to 100 KHz. 
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Figure 8. Enlarged photograph of the chip. 
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Figure 9. Magnitude response of the 40 Hz design with 10 kHz clock. 
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Figure 10. Phase response of the 40 Hz design with 10 kHz clock. 
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Figure 1 1 . Magnitude response of the 2.5 Hz design with 160 kHz clock. 
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Figure 12. Phase response of the 2.5 Hz design with 160 kHz clock. 
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